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Abstract We present a new experimental system (the "atom-optics billiard" ) and 
demonstrate chaotic and regular dynamics of cold, optically trapped 
atoms. We show that the softness of the walls and additional optical 
potentials can be used to manipulate the structure of phase space. 
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1. Introduction 

The possibility of using lasers to manipulate, cool and trap neutral 
atoms has revolutionized many areas of atomic physics, and opened new 
research fields. In particular, cold atomic samples trapped in far-detuned 
optical traps are an ideal starting point for further experimental work, 
since the atoms can be confined to a nearly perturbation-free environ- 
ment, Doppler shifts are almost suppressed, and long interaction times 
are possible (Grimm et al., 2000). 

Cold atoms have been used before in non-linear dynamics experi- 
ments. First, by confining atoms to the minima of a time-dependent 
amplitude (or phase) periodic optical potential, effects such as dynamical 
localization were observed for ultra-cold sodium atoms (Raizen, 1999). 
More recently, a similar system was used to observe chaos- assisted tun- 
nelling (Steck et al., 2001; Hensinger et al., 2001). 

The "billiard" (a particle moving in a closed region of coordinate 
space and scattering elastically from the surrounding boundary) is a 
well known paradigm of non-linear dynamics and much theoretical work 
has been done about its properties. An interesting experimental mani- 
festation of classical billiard dynamics is the observation of conductance 
fluctuations in semiconductor microstructures in the coherent ballistic 
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regime (Jalabert et al., 1990; Marcus et al., 1992; Baranger and Wester- 
velt, 1998). The structure of this fluctuations is affected by the classical 
dynamics of the electrons. Chaotic dynamics were also exploited to con- 
struct lasers with a very high power directional emission (Gmachl et al., 
1998). 

We present a new experimental system, the "atom-optics billiard" 
in which a rapidly scanning and tightly focused laser beam creates a 
time-averaged quasi-static potential for ultra-cold Rb atoms (Friedman 
et al., 2000). By controlling the deflection angles of the laser beam, we 
create various billiard shapes. As opposed to textbook billiards and other 
experimental realizations, the particles in our system have an internal 
quantum structure. The internal degrees of freedom are coupled to the 
external ones, and spectroscopy can serve as a sensitive probe for the 
dynamics (Andersen et al., 2003). 

In this lectures, we review several classical non-linear dynamics ex- 
periments in atom-optics billiards. The dynamics of the atoms confined 
by such billiards is studied by measuring the decay in the number of 
confined atoms through a small hole on the boundary. We demonstrate 
the existence of regular and chaotic motion, and study the effects of scat- 
tering by impurities on the dynamics (Friedman et al., 2001). Next, we 
study the changes in the structure of phase-space induced in a chaotic 
billiard when the billiard's boundaries are soft and find that in a chaotic 
Bunimovich stadium soft walls induce islands of stability in phase-space 
and greatly affect the system as a whole (Kaplan et al., 2001; Kaplan 
et al., 2004). Finally, we find that introducing additional potentials 
can stabilize an otherwise chaotic billiard, or alternatively destabilize an 
otherwise regular one, by inducing curvature in the particles trajectories 
between bounces from the billiards walls (Andersen et al., 2002). 

2. Optical Dipole Traps 

The interaction of a neutral atom and a light field is governed by 
the dipole interaction, and is usually separated into two terms which 
correspond to a reactive force and a dissipative force. When an atom is 
exposed to light, the electric field component E induces a dipole moment 
d in the atom, oscillating at the driving light frequency. The amplitude 
of the dipole moment is related to the amplitude of the field by d = aE, 
where a is the atomic complex polarizability, which is a function of the 
driving frequency. The interaction of the induced potential with the 
driving field gives rise to the potential: 
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U dip = -\ (d- E) = -2-Re( a )m, (1) 

2 

where / = 2eoc|.E| is the light intensity. The reactive "dipole" force is 
a conservative one, and is equal to the gradient of the above potential. 
The dissipative force is related to the power the oscillating dipole absorbs 
from the field, which is given by: 

P ohs = (d-E) = — Im(a)/(f). (2) 
\ / e c 

In a quantum picture, the dipole force results from absorption of a 
photon from the field followed by stimulated emission of a photon into 
a different mode of the laser field. The momentum transfer is the vector 
difference between the momenta of the absorbed and emitted photons. 
The dissipative component has its origin in cycles of absorption of pho- 
tons, followed by spontaneous emission, in a random direction. Using 
Eq. 2 we can write then an equation for the rate of spontaneous photon 
scattering: 

The atomic polarizability can be calculated by using the solutions of the 
optical Bloch equations, while translational degrees of freedom are taken 
into account (Cohen- Tannoudji, 1990). Using this result for a two- level 
atom, introducing the "Rotating Wave Approximation" , and assuming 
5 3> 7, the dipole potential and spontaneous photon scattering rate can 
be written as: 

3ttc 2 7 



= %4i' (r> (4) 



where fkJo is the energy separation of the atomic levels, 5 = to — too the 
laser detuning and 7 the natural linewidth of the transition. 

The dissipative part of the atom-light interaction is used for laser 
cooling of atoms (Cohen- Tannoudji, 1990), a pre-requisite for optical 
trapping. However, spontaneous photon scattering is in general detri- 
mental to trapped atoms, mainly since it can induce heating and loss. 
Comparing the expressions for the dipole force and scattering rate, under 
the above approximations, yields the relationship 
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which indicates that a trap with an arbitrarily small scattering rate can 
be achieved by increasing the detuning while maintaining the ratio 1/5. 

Equations 4 and 5 indicate that if the laser frequency is smaller than 
the resonance frequency, i.e. 5 < ( "red-detuning" ) the dipole potential 
is negative and the atoms are attracted by the light field. The minima 
of the potential is found then at the position of maximum intensity. In 
the case 5 > ( "blue detuning" ) the minima of the potential is located 
at the minima of the light intensity. 

Trapping atoms with optical dipole potentials was first proposed by 
Letokhov (Letokhov, 1968) and Ashkin (Ashkin, 1970) . Chu and cowork- 
ers (Chu et al., 1986), were the first to realize such a trap, trapping about 
500 atoms for several seconds using a tightly focused red-detuned beam. 
Later, a far-of-resonant-trap for Rb atoms was demonstrated (Miller 
et al., 1993), with a detuning of up to 65 nm, i.e. 5 > 5 x 10 6 7. In this 
case, the potential is nearly conservative and spontaneous scattering of 
photons is greatly reduced. A comprehensive review of the different 
schemes and applications of such optical dipole traps is presented in 
(Grimm et al., 2000). In the limiting case where the frequency of the 
trapping light is much smaller then the atomic resonance, trapping is 
still possible, practically with no photon scattering (Takekoshi et al., 
1995). Such a quasi-electrostatic trap, formed by two crossed CO2 laser 
beams, was used to create a Bose-Einstein condensate without the use 
of magnetic traps (Barrett et al., 2001). 

Apart from using far-off-resonance lasers, the interaction between the 
light field and the atoms can be reduced by the use of blue-detuned 
traps, in which atoms are confined mostly in the dark. In the first dark 
optical trap cold sodium atoms where trapped using two elliptical light 
sheets which where intersected and formed a " V" cross-section (Davidson 
et al., 1995). Confinement was provided by gravity in the vertical direc- 
tion and by the beams' divergence in the longitudinal direction. Many 
other configurations where proposed since (Friedman et al., 2002). Of 
special importance are traps formed with a single laser beam, which are 
particularly simple to align, hence it is easier to optimize their proper- 
ties. The simplicity of these traps also enables easy manipulation and 
dynamical control of the trapping potential, its size and its shape. 

3. Experimental Realization of Atom-Optic 
Billiards 

Friedman and coworkers (Friedman et al., 2000) developed a new 
dark optical trap for neutral atoms: the ROtating Beam Optical Trap 
(ROBOT). In the ROBOT a repulsive optical potential is formed by a 
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Figure 1. Examples of atom-optic billiards of different shapes. Shown are CCD 
camera images at the rotating beam focus. From left to right: A circular billiard, an 
elliptical billiard, a tilted Bunimovich stadium with a hole, and a Sinai billiard with 
a hole. 



tightly focused blue-detuned laser beam which is a rapidly scanned using 
two perpendicular acousto-optic scanners (AOS's). The instantaneous 
potential is given by the dipole potential of the laser Gaussian beam: 



U (x,y,t) 



k 7) exp 



-2^x-xo {t)f + {y-yo (i)) 2 ) /w 2 



(7) 



where (xq (t) , yo (t)) is the curve along which the center of the Gaussian 
laser beam scans, wo is the laser beam waist, P is the total laser power, 
and k = i^tj- When the x and y scanners perform a sinusoidal and 

cosinusoidal scan, respectively, the time-averaged intensity in the focal 
plane is a ring with a Gaussian cross-section. The radial potential barrier 
is then given by 



U = k (-) 1/2 — — , (8) 

\TT J ZTtrWo 

where r is the scanning radius. 

Several measurements were performed to prove that the potential can 
be regarded as a time averaged potential (e.g. measurements of the 
oscillation frequencies of atoms in the trap, and the stability as a function 
of the scanning frequency). 

We used the ROBOT to create what we denoted "Atom-Optics Bil- 
liards" . By controlling the deflection angles of both AOS's using arbi- 
trary waveform generators, arbitrary scans are performed independently 
in the x and y direction, thus creating an arbitrary 2D scan of the beam 
in the focal plane. We create various billiard shapes, such as an ellipse, 



1 Optical billiards for atoms were developed also in the Raizen group, simultaneously with 
our work (See (Milner et al., 2001)) 
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a circle and a tilted Bunimovich stadium that confine the atoms in the 
transverse direction (See Fig. 1), and probe their dynamical properties 
by measuring their decay curve through a hole in their boundary. 

The departure of our system from the mathematical definition of a 
billiard due to interactions, softness of the walls, additional forces acting 
on the atoms, etc., is probably what makes it interesting, and we will 
devote most of these lectures to these effects. But first, we would like to 
explain why, in spite of all the latter, we can treat our optical traps as 
textbook billiards. 

Three dimensional structure: The effects we study are two di- 
mensional effects. Fortunately, a large difference exists in time scales 
for the dynamics in the radial and longitudinal dimensions and enables 
us to treat the billiards as 2D billiards, provided than our experiments 
are done in a short enough time. Alternatively, to better approximate 
a true two-dimensional system, we confine the atoms by a stationary 
blue-detuned standing wave along the optical axis. Here, the atoms are 
tightly confined near the node-planes of the standing wave, but move 
freely in these planes, forming "pancake" shaped traps separated by ~ 
400 nm (half the wavelength of the standing wave laser). 

Velocity distribution: The trap is loaded from a thermal cloud 
of atoms, and the trapped atoms have a nearly thermal distribution of 
velocities with a root-mean-square (RMS) width of ~ llv reco u, and a 
typical velocity similar to this width (v reco u is the velocity acquired by 
an atoms that recoils after absorbing a single photon, 6 mm/s for the D2 
line in 85 Rb atoms). To better approximate a mono-energetic case, where 
all atoms have identical velocity, we expose the atoms, after loading into 
the billiard, to a short (1.5 us) pulse of a strong, on-resonance, "pushing" 
beam perpendicular to the billiard beam and at 45° to the vertical axis. 
Following this pushing beam, the center of the velocity distribution is 
shifted to 20v reco n, while the RMS width grows to l2v reco n. After an 
additional 50 ms of collisions with the billiard's boundaries, the direction 
of the transverse velocity distribution is completely randomized by the 
time the hole is opened. 

Gravitation: Typically, the mean gravitational energy in the traps 
is \mgh < 70E reco n, where h is the vertical size of the billiard and 
Erecoil = k mv recoii * s * ne recou energy. After the pushing beam, the 
atom's kinetic energy is ~ 400E recoi i and therefore gravity can be con- 
sidered a small perturbation. In some cases, this perturbation can be of 
great importance (as seen in 10.2) but in general it can be neglected. 

Collisions: In the range of atomic densities realized in the billiard, 
the mean collision time between atoms is much longer than the experi- 
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(a) (b) (c) 

Figure 2. Numerical simulation of two dimensional trajectories of Rb atoms in 
different atom-optics billiards, (a) Elliptical billiard: Only one of the two existing 
types of trajectories is shown, for which atoms are confined by a hyperbolic caustic, 
and thus excluded from a certain part of the boundary, (b) Circular billiard: Nearly 
periodic trajectories demand an increasingly long time to sample all regions on the 
boundary, (c) "Tilted" stadium: Every atomic trajectory would reach a certain region 
in the boundary with a comparable time scale. 



merit time, hence the motion of the atoms between reflections from the 
walls can be regarded as strictly ballistic. 

Spontaneous Photon Scattering: Scattering events involve a change 
in momentum and can cause a diffusion in phase space which is destruc- 
tive to any dynamical effect. To minimize photon scattering and thus 
ensure ballistic motion of the atoms and elastic reflections from the bil- 
liard potential walls, we use far-detuned laser beams to form the billiard 
(5 = 0.5 — 1.5 nm > 4 • 10 4 7) and standing wave (5 = 4 nm > 3 • 10 5 7). 

We performed also numerical simulations in order to check the above 
assumption of billiard-like dynamics. The simulations include all our 
experimental parameters: a Gaussian beam scanning at a finite rate, 
a thermal ensemble, 3D structure of the beam and atomic cloud, the 
loading process. Shown in Fig. 2 are typical trajectories of a single 
atom (the simulation shown here is, for illustrative purposes, a 2D one). 
The motion of the atom shows the same features of the motion in an 
ideal billiard: Periodic or quasi-periodic for the elliptical (Fig. 2a) and 
circular (Fig. 2b) billiard and chaotic for the Bunimovich "stadium" (Fig. 
2c). 

Our experimental procedure is described in (Friedman et al., 2001). 
Briefly, ~ 3 x 10 8 atoms are loaded, during 700 msec, from the vapor 
cell into a magneto-optical trap (MOT) (Raab et al., 1987). Next, after 
47 ms of "weak" and "temporally dark" MOT (Adams et al., 1995), 
and 3 ms of polarization gradient cooling (PGC) (Dalibard and Cohen- 
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Tannouji, 1989), during which the trapping beam is on, we end up with 
0.5 — lxlO 6 atoms in the billiard with a temperature of ~ 10/iK and 
a peak density of ~5xl0 10 cm -3 . After all the laser beams are shut 
off (except for the rotating beam), a hole is opened in the billiard po- 
tential, by switching off one of the AOS's for ~ 1 fis every scan cycle, 
synchronously with the scan. The number of atoms remaining in the 
trap is measured using fluorescence detection, with a 100/us pulse reso- 
nant with \5Si/2,F = 3) — > |5P 3 / 2 ,-F = 4). The ratio of the number of 
trapped atoms with and without the hole, as a function of time, is the 
main data of our experiments. 

4. Integrable and Chaotic Systems 

The billiard is a conservative Hamiltonian system. We discuss here 
some general features of this class of systems. The equations of motion 
for an ^-dimensional Hamiltonian system are: 



Qk 



d Pk Pk ~ dq k ' 



with k = 1 . . . N. If there are N independent, "well-behaved" inte- 
grals Fi(q,p), that are constant along each trajectory, then the system 
is called Integrable. For a conservative system, one of this constants is 
the energy. The trajectory of an integrable system in phase-space lies 
in the intersection of Fi(q,p) = fi where fi are constants. Hence, the 
existence of N integrals of motion implies that the trajectory is confined 
to an iV-dimensonal manifold, A4, in the 2iV-dimensonal phase-space. 
This iV-dimensional manifold does not have necessarily a simple struc- 
ture, but if the integrals of motion are "in involution" , i.e. their Poisson 

brackets with each other vanish i^j =0, then M. has the topology 

of an N-dimensional torus (Schuster, 1984; Gutzwiller, 1990). 

A natural way to describe the dynamics is by using a set of coordi- 
nates that define (i) on which torus the trajectory takes place, and (ii) 
a coordinate on the torus. The standard way to do that is by using 
a canonical transformation to "Action and Angle" coordinates, (1,0), 
which are defined as those variables that transform TC(q,p) into T~C(I), 
i.e. a Hamiltonian which does not depend on one of the new variables. 
Such a transformation always exists for an integrable system (Gutzwiller, 
1990). The equations of motion will be then: 

h = - d -0 = o k = ™& = . k (I), (10) 
where u) is a constant. The equations are then trivially integrated giving: 
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I(t) = 1(0) 9(t) =cot + 9(0) . (11) 

I are the actions of the torus, and 9 the angles on the torus. Hence, 

uj = 9 are the "frequencies" on the torus. Close ("periodic") orbits occur 
only if the frequencies have a rational ratio, i.e. — = — with m, n £ N, 
for any i,j = 1...N. For irrational frequency ratios the orbit never 
repeats itself but approaches any point on the manifold infinitesimally 
close in the course of time (we will call this trajectories "quasi-periodic" ) . 

Irregular or non-integrable trajectories fill a subspace with grater di- 
mensionality, possibly with the same dimensionality as the whole phase 
space. Some of this irregular regions of phase space are "stochastic", 
meaning that, although they are deterministic, they display an extremely 
sensitive dependence on the initial conditions and are unpredictable in 
the long time. The field of nonlinear dynamics is rich in ways to char- 
acterize a system which exhibits stochastic dynamics. For example, an 
"ergodic" system is one in which time averages (over a single trajectory) 
are equal to ensemble averages. This implies that every single trajec- 
tory "samples" the entire phase space. The extreme case of stochasticity 
are K-systems or systems exhibiting "strong chaos", which are defined 
as systems in which initially close orbits separate exponentially. Terms 
such as Mixing and positive K-entropy are also commonly used. We 
will not be concerned with this fine-tuned hierarchy and will denote our 
irregular systems, in a loose way of speaking, "chaotic". 

It should be noted, that circular, rectangular and elliptical billiards 
(which exhibit integrable motion), and Sinai and Bunimovich billiards 
(in which the dynamics is completely chaotic), are the exemptions in 
Nature. Most systems have a "mixed" phase space in which regular and 
chaotic regions coexist. 

5. Poincare Surface of Section 

Let us restrict ourselves to a 2-dimensional system. Phase space is 
four-dimensional, and the TL = E surface has three dimensions. A con- 
venient way to visualize the dynamics is by using a "Poincare Surface 
of Section" obtained by recording the successive intersections of a tra- 
jectory with an arbitrary plane of section. For example, the x surface 
of section S x is the intersection of the energy surface with y = and 
has (p x , x) as coordinates (see Fig. 3a) . The value of p y is determined 
by the condition TC(x,y = 0,p x ,p y ) = E, up to a sign. If we define it 
as positive, we see that specifying the position of a system on S x com- 
pletely specifies its state. An initial point Xq = (xo,p XQ ) will cross S x 
again at X\ = (xi,p Xl ), X2 = (x2,p X2 ) and so on. The resulting map 
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Figure 3. Two ways to construct a Poincare surface of section, (a) Phase space 
is represented using x (the points at which the orbit crosses the y — line), and 
p x = v cos#. (b) Bouncing Map: the usual coordinates are the arc length, s, and the 
tangential momentum p — cos a. 



Xi = T(JQ_i) of S x onto itself (the Poincare map) is a simplified picture 
of the dynamics. 

Another common way to construct a Poincare section for a billiard (see 
Fig. 3b), is by looking at its bouncing map, which specifies the evolution 
of position and momentum from one collision with the boundary to the 
next one. The position on the boundary can be parameterized either by 
the arc length s or the direction of the tangent, ip, and the momentum 
can be specified using the particle's direction with respect to the tangent, 
a, or the tangential momentum, i.e. p = cos a (the coordinates (s,p) 
are called Birkhoff coordinates). 

Repeated iterations of the Poincare map reveal whether or not the mo- 
tion is integrable. In the case of integrable motion, the system explores 
only a 2D torus, and its intersection with S is a smooth, closed curve, 
which becomes apparent after enough iterations. If the orbit closes (i.e. 
the motion is periodic) then X n = Xq for some n, depending of the ratio 
The torus in this case is full of such curves, and the curve will appear 
when the iterations are made for many different starting points Xq. In 
the more common case of irrational ^ the curves will be generated by 
iterating any initial point Xq a large enough number of iterations. For 
non- integrable motion tori do not exist, and the system explores a three 
dimensional region. Its crossings with S will not cover a curve, but a 
two dimensional region of S. 

6. Decay through a hole 

One way to characterize the dynamics in a billiard is by looking at 
the escape probability of a particle from it through a small hole in its 
boundary or, alternatively, at the decay in the number of particles in it 
as a function of time. This decay rate is strongly related to the decay 
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in the system correlations, which has been a topic of numerous studies 
(See for example (Friedman and Martin, 1984)). A different functional 
dependence is expected for chaotic and regular billiards. 

Assume that a billiard with area A contains an ensemble of particles 
with momentum in the range \p,p+Sp], and that there is a small hole, of 
length /, in the billiard's boundary. It can be shown (Bauer and Bertsch, 
1990) that if the billiard is completely ergodic, the number of particles 
in it, as a function of time, obeys 

7T A 

N(t) = N(t)/r c , r c = — . (12) 

Hence, the decay from a chaotic and ergodic billiard exhibits a universal 
(i.e. similar for all chaotic billiards) exponential behavior, with a decay 
time constant r c . The decay from an integrable billiard does not have a 
universal functional dependence, and differs for different billiard shapes. 
Nevertheless, power-law decays are very common. 

It is important to note already at this point that even for "completely 
chaotic" systems there can be deviations, in the long time tail of the de- 
cay curve, from the above exponentiallity. For the Bunimovich billiard, 
for example, an algebraic tail has been observed, when the hole in not 
infinitesimalry small (Legrand and Sornette, 1990; Alt et al., 1996). This 
tail is induced by the marginally stable "bouncing ball" orbits, where 
a particle bouncing back and forth between the opposite straight lines. 
An orbit near one of these closed orbits, will perform a "zig-zag" path 
for many bounces (a "resonance") before getting lost in the chaos and 
the particles can "stick" for a long time to the vicinity of the boundary 
of the marginally stable region. The phenomena of "stickiness" will play 
an important role in 9. Note that in our experiments, we choose to used 
a "tilted" Bunimovich stadium, in which the two semi-circles have differ- 
ent size and the straight lines are not parallel, in order to preclude this 
strong position-dependent correlations between bounces (Vivaldi et al., 
1983). 

7. Chaotic and Integrable Dynamics 

The radius of curvature of a circular billiard is constant along the 
boundary, and the orbit of a particle in it is composed of a succession 
of chords making angles a with the boundary's tangent (Berry, 1981). 
The tangential momentum, p = cos a, is then conserved and the system 
is integrable. The well-known Bunimovich stadium is composed of two 
semicircular arcs with radius R joined by straight lines with length L. 
For any < L < oo the system is chaotic (Bunimovich, 1974). As 
already mentioned above, we use a "tilted" Bunimovich stadium, in 
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Figure 4- Decay of atoms from circular and stadium billiards: The decay from the 
stadium billiard (•) shows a nearly pure exponential decay. For the circle (♦) the 
decay curve flattens, indicating the existence of nearly stable trajectories. The full 
lines represent numerical simulations, including all the experimental parameters, and 
no fitting parameters. The dashed line represents exp(-t/r c ), where r c is the escape 
time calculated for the experimental parameters. 



which the two semi-circles have different size and the straight lines are 
not parallel. 

Comparing the atomic trajectories for the circular billiard (Fig. 2b) 
and the tilted-stadium billiard (Fig. 2c) illustrates what can be de- 
noted a microscopic effect in phase space. Neither of these shapes has 
a macroscopic stable region in phase space for a hole at any point on 
the boundary. Nevertheless, as explained in 6, differences in the decay 
rate are expected. For the circular billiard nearly periodic trajectories 
exist that require an arbitrarily long time to sample all regions on the 
boundary (the exactly periodic trajectories that are completely stable 
have only a zero measure, and hence can be neglected). This yields 
many time scales for the decay through a small hole on the boundary, 
and results in an algebraic decay (Bauer and Bertsch, 1990). For the 
tilted-stadium billiard, phase space is chaotic resulting in a pure expo- 
nential decay, with a characteristic decay time r c (Bauer and Bertsch, 
1990). 

The experimental decay curves for the circular billiard and the tilted- 
stadium with identical potential height, area and hole length, and with 
the hole located at the bottom, are shown in Fig. 4. The area of the 
circular billiard and the stadium were set to be equal, with a precision of 
~ 10%. The initial decay is nearly identical for the two shapes. This is 
expected since it takes a certain time for the atomic ensemble to "feel" 
the shape of the billiard. However, for longer times the decay curve for 
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the circular billiard flattens, indicating the existence of nearly stable tra- 
jectories, whereas that of the stadium remains a nearly pure exponential. 
Also shown in Fig. 4 are the results of full numerical simulations that 
contain no fitting parameters. The simulations include the measured 
three dimensional atomic and laser-beam distributions, atomic velocity 
spread, laser beam scanning, and gravity. As seen, the simulations fit 
the data quite well. For long times the circular billiard simulations pre- 
dict greater stability than the experimental data. A possible explanation 
for this difference are imperfections of the billiard shape. Including im- 
perfections in the beam-shape in the simulations yielded similar effects 
of reducing the long-time stability of the circular billiard. The devi- 
ation of the stadium decay curve from the simulation results, can be 
explained by the existence of a small amount of atoms trapped in areas 
far from the focus, where the real intensity distribution deviates from 
the simulated one to a greater extent. Finally, Fig. 4 also shows a pure 
exponential decay curve with a time constant r c = 4.9 ms, calculated 
with the measured billiard parameters. As seen, it closely resembles the 
full simulation results, indicating that for this configuration, the stadium 
behaves nearly like an ideal chaotic billiard (two-dimensions, no gravity, 
mono-energetic, zero wall thickness). Note, that deviations from a pure 
exponential decay, representing correlations due to the finite hole size, 
are expected at long times (Alt et al., 1996; Vivaldi et al., 1983) but 
occur below the noise level of our experiment. 

8. Elliptical Billiard and the Effect of Scattering 
by Impurities 

The motion in an elliptical billiard is integrable (the integral of motion 
is the product of the angular momenta about the two foci (Berry, 1981)). 
It can also be shown (Koiller et al., 1996) that if one segment of the 
trajectory crosses the line connecting the two foci, then every segment 
will cross it, and hence there are two families of trajectories: First, 
"external" trajectories are those that never cross the line between foci. 
They bounce all round the billiard, and are confined outside elliptical 
caustics, smaller than the billiard itself but with the same focal points. 
The second family, "internal" trajectories, are those that cross the line 
connecting the two focii. They explore a restricted part of the boundary 
and are confined by hyperbolic caustics, again with the same focal points. 
Phase space is then macroscopically divided into two separate regions 
corresponding to the types of orbits described above (Koiller et al., 1996). 

The decay in a system with such a non-uniform phase space, is very 
sensitive to the position of the hole. If it is located at the short side of the 
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Figure 5. Decay of the number of atoms from an elliptical billiard: The full symbols 
denote the unperturbed case, in which the surviving fraction for the ellipse with the 
hole on the long side (■) decays much faster than for the hole on the short side (•). 
The insets show CCD camera images of the billiards' cross-sections at the beam focus. 
The empty symbols show the case in which a 10 /is velocity randomizing molasses 
pulse is applied every 3 ms (see text). 



ellipse, particles in internal trajectories will remain confined and never 
reach the hole. Alternatively, all trajectories, excluding a zero-measure 
amount, reach the vicinity of a hole on the long side of the ellipse and 
hence the number of confined particles decays indefinitely. 

Figure 5 shows the measured survival probability for the elliptical bil- 
liard with a 60 /im hole on the long side and on the short side. These 
measurements are taken without a pushing beam. To minimize the effect 
of gravity, the ellipse is rotated such that the hole's direction is perpen- 
dicular to gravity for both cases (see inset). The results reveal that the 
initial decay rate (for the first few points) is identical for both cases, as 
expected, confirming the experimental accuracy of our shapes and holes 
sizes. However, at longer times, the survival probability for the hole on 
the short side becomes much higher than for the hole on the long side, 
as expected from the discussion above. In general, the stability due to 
such a macroscopic phase-space separation is very robust and remains 
nearly unchanged with or without the standing wave, for a large variety 
of hole sizes, for different orientations of the ellipse relative to gravity, 
and with small distortions of the ellipse's potential shape. 

Next, we introduce a controlled amount of scattering to the atomic 
motion. We expose the confined atoms to a series of 10 (is pulses of PGC 
(using the six MOT beams) every 3 ms. During each pulse, each atom 
scatters ~30 photons, and hence its direction of motion is completely 
randomized, whereas the total velocity distribution remains statistically 
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Figure 6. The surviving fraction of atoms 150 ms after the hole was opened on the 
long (□) and short (o) side of the elliptical billiard of Fig. 5, as a function of the 
repetition rate of the velocity randomizing molasses pulses. At a rate comparable to 
the escape rate, r c ~ , the stability of the trajectories was significantly reduced. 



unchanged. Moreover, the atoms barely move during each short pulse, 
and maintain ballistic motion between the pulses. Thus, the effect of the 
pulses resemble the effect of random scattering from fixed points (impu- 
rities) or of short range binary atomic collisions. The measured decay 
curves for this case are also shown in Fig. 6, for the two hole positions 
of the ellipse. As seen, for the hole on the long side the randomizing 
molasses pulses cause little change. However, for the hole on the short 
side, a complete destruction of the stability occurs, and the decay curves 
for the two hole positions approximately coincide. 

This effect is further illustrated in Fig. 6, where the measured sur- 
vival probability (150 ms after the hole was open) as a function of the 
repetition rate of the velocity-randomizing pulses, is shown for the two 
hole positions in the ellipse. As seen, the pulse rate required to signifi- 
cantly reduce the stability is approximately one pulse per r c , calculated 
as 12 ms here, when averaging over the thermal velocity distribution of 
the atoms. 

9. Billiards with Soft Walls 

Kolmogorov-Arnold-Moser (KAM) theory provides a framework for 
the understanding of an integrable Hamiltonian system which is subject 
to a perturbation (Gutzwiller, 1990). As the strength of the pertur- 
bation is increased areas of phase-space become chaotic and eventually 
only islands of stability, surrounded by a chaotic "sea" , survive. The op- 
posite question is of interest as well: What will happen to a completely 
chaotic system when a perturbation is applied? will it remain chaotic or 
will islands of stability appear?. For the billiard system, one such per- 
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turbation is making the walls "soft" , as opposed to the infinitely steep 
potential of the ideal case. This kind of perturbation is also interesting 
in the context of understanding the origins of statistical mechanics, for 
which the billiard problem is a widely used paradigm (Zaslavsky, 1999). 
The Sinai billiard (Sinai, 1970) is mathematically analogous to the mo- 
tion of two disks on a two-dimensional torus, an approximation for gas 
molecules in a chamber, and has been proved to be ergodic. However, if 
the potential between the two disks is smooth, as is the actual potential 
between gas molecules, it was shown that elliptic periodic orbits exist, 
hence the system is not ergodic (Donnay, 1999). 

From the practical point of view, physically realizable potentials are 
inherently soft, and the softness of the potential may result in a mixed 
phase-space with a hierarchical structure of islands. This structure 
greatly affects the transport properties of the system (e.g. induces non- 
exponential decay of correlations), since trajectories from the chaotic 
part of phase-space are trapped for long times near the boundary be- 
tween regular and chaotic motion (Zaslavsky, 1999). As an important ex- 
ample, these considerations were studied in the context of ballistic nano- 
structures, and found to be the cause for the fractal nature of magneto- 
conductance fluctuations in quantum dots (Ketzmerick, 1996; Sachrajda 
et al., 1998). However, in these systems the wall softness is often ac- 
companied by non-ideal effects (such as scattering from impurities) and 
hence its role is still controversial. 

For a certain kind of dispersing billiards it was theoretically proven 
that when the wall becomes soft, an island appears near a singular, tan- 
gent trajectory (Turaev and Rom-Kedar, 1998; Rom-Kedar and Turaev, 
1999). More recently, and partly inspired by our numerical results, the 
appearance of an additional island, this time near a corner, was also 
shown (Turaev and Rom-Kedar, 2003). 

The first billiard that we study is a tilted Bunimovich stadium (see 
insets in Fig. 7), composed of two semicircles of different radii (64 ^m 
and 31 /mi), connected by two non-parallel straight lines (192 ^m long). 
When the potential of the wall becomes softer, a stability region appears 
around the singular trajectory which connects the points where the big 
semicircle joins the straight lines. 

We use the laser beam 1/e 2 radius, wo, to set the softness of the 
billiard's walls (See Eq. 7), and experimentally control it by the use of a 
telescope with a variable magnification, located prior to the AOS's such 
that wq could be changed without affecting the billiard's size and shape. 
As seen from Eq. 8, the time averaged potential is kept constant for 
different values of wq by changing the laser power. 
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Figure 7. Experimental results for the decay of cold atoms from a tilted-stadium 
shaped atom-optics billiard, with two different values for the softness parameter: 
wq — 14.5pim (o), and wq = 24/im (+), and for two different hole positions, (a) 
The hole is located inside the big semicircle. The smoothening of the potential wall 
causes a growth in stability, and a slowing down in the decay curve, (b) The hole 
includes the singular point, no effect for the change in wo is seen. Also shown are 
results of full numerical simulations, with the experimental parameters (see text) and 
no fitting parameters. The dashed line is e~ t ^ Tc , the decay curve for an ideal (hard- 
wall) billiard. The insets show measured cross-sections of the (averaged) intensity of 
the laser creating the soft-wall billiards, in the beam's focal plane. The size of the 
images is 300 x 300/im. 



In Fig. 7, experimental results for the decay from a tilted stadium 
with two different values of the softness parameter (wo = 14.5 fim and 
wo = 24 ^m), are presented. It can be seen that when the hole is 
located entirely inside the big semicircle (Fig. 7a), the soft wall causes 
an increased stability, and a slowing down in the decay curve. When the 
hole includes the singular point where the semicircle meets the straight 
line (Fig. 7b), no effect for the change in wq is seen. We show below 
that these results can be explained by the formation of a stable island 
around the singular trajectory, and a sticky region around it. Figure 
7 includes also the results of numerical simulations, which include the 
three dimensional atomic and laser-beam distributions, atomic velocity 
spread, laser beam scanning and gravity, and no fitting parameters. As 
can be seen, there is a very good agreement between the simulated and 
measured decay curves. Similar decay measurements and simulations for 
a circular atom-optics billiard show no dependence on wq in the range 
14.5 — 24 //m, and no dependence on the hole position. 

To understand these observations, it is useful to look on how the 
phase-space of the system changes with changing wo. In Fig. 8, results 
of numerical simulations for classical trajectories of Rb atoms inside the 
tilted-stadium billiard are shown. For clarity, we assume a monoener- 
getic ensemble (with v = 20v reco n), a two-dimensional system, and no 
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Figure 8. Poincare surface of section for monoenergetic atoms confined in a tilted- 
stadium atom-optics billiard with parameters specified in text, and three different 
values of the softness parameter wo- (a) wo = 18/im. A small elliptic island appears 
close to the trajectory which connects the two singular points. Upper inset: a tra- 
jectory in the island. Lower inset: a typical chaotic trajectory in this billiard, (b) 
wo = 27/nm. Three additional islands appear around the central one, and correspond 
to the periodic trajectory shown in the upper inset. Around these islands there is a 
large area of "stickiness", where the trajectories spend a long time (see trajectory in 
lower inset), (c) wo = 30/im. The three previous islands merge into one large elliptic 
island (see trajectory in upper inset), with some stickiness around it (lower inset). 



gravity. The dimensions of the billiard are equal to the experimental 
ones. Phase-space information is presented using a Poincare surface of 
section, showing v x versus x at every trajectory intersection with the 
billiard's symmetry axis (y = 0), provided that v y > 0. 

In Fig. 8(a), the wq = 18 /jm case is shown. A small elliptic island 
appears around the trajectory which connects the two singular points, 
as can be seen also from the upper inset, which shows a trajectory in 
the island. In the lower inset, a typical chaotic trajectory is shown. 
For wo < 12 jum, no islands with area larger than 10~ 4 of the total 
phase-space (the resolution of the simulations) are observed. In general, 
the island size increases with the increase of wq, as can be seen from 
Fig. 8(b), (c) which correspond to wo = 27,30 /im, respectively. For 
wq = 27 /im, three additional islands appear around the central one, 
and correspond to the periodic trajectory shown in the upper inset. 
Around these islands there is a large area of "stickiness" 2 , where the 
trajectories spend a long time. Such a "sticky" trajectory is presented 
in the lower inset of Fig. 8(b). The exact structure of the island and its 
vicinity depends on the softness parameter wo in a sensitive way, as can 
be seen from the wq = 30 /im case (Fig. 8(c)), where the three previous 



2 We define a "sticky" trajectory as one which spends inside a box surrounding the island 
(shown in Fig. 2) a time which is more than x3 longer than expected for a random trajectory. 
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Figure 9. Fraction of remaining atoms at 5t c after the hole is open, as a function 
of the softness parameter (wo), for a hole in the big semicircle. (+): experimental 
results, (o): numerical simulation. The dashed line is e~ 5 , the expected value for 
small wo- Also shown are values for the island (a) and island + stickiness (•) sizes as 
a fraction of the phase-space area, calculated from the two-dimensional phase-space 
simulations. 

islands merge with the central one into one big elliptic island, with some 
stickiness around it. 

In Fig. 9, the measured fraction of remaining atoms at 5r c (= 42.5 
ms) after a hole is opened in the big semicircle is plotted as a function 
of the softness parameter, wq, together with the results of numerical 
simulations. A very good agreement exists between the decay simula- 
tions and the measured data, and both converge to the expected value 
of e~ 5 for small wq. To intuitively understand the origin of the increased 
stability, the results of numerical simulations, showing the relative area 
in phase space of the island and the combined area of the island and 
the "sticky" regions, are also presented. The simulations reveal that the 
total area size of the island and sticky regions grows monotonically with 
wq, in a similar way to the decay results. The island size itself is much 
smaller, and its size has a non-monotonic dependence on wq. These facts 
suggest that the remaining atoms in the experiment, at 5r c , are mainly 
due to stickiness and demonstrate the important effect of the stickiness 
on the dynamics of the billiard. 



If particles in the billiard move in straight lines between bounces from 
the wall, then the movement between them is neutral, and the dynamics 
is determined solely by the shape of the boundary. Exponential sep- 
aration of close trajectories (chaotic dynamics) can only occur if the 
boundary introduces sufficient instability. But if particles move along 
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curved trajectories between bounces from the wall the dynamics cru- 
cially depends on the specific properties of this motion (we call these 
billiards "curved-trajectory billiards"). Curved trajectory billiards have 
been found relevant for micro devices (Takagaki and Ploog, 2000), and 
much effort have been put into studying them. In these studies the curva- 
ture of the trajectory arises from a magnetic field (Tasnadi, 1997; Robnik 
and Berry, 1985; Gutkin, 2001; Dullin, 1998) or from the curvature of 
the surface on which the billiard is made (Gutkin et al., 1999). Billiards 
where gravity provides the confining force along its axis have been stud- 
ied experimentally, theoretically and numerically (Lehtihet and Miller, 
1986; Wallis et al., 1992; Hughes et al., 2001; Milner et al., 2001). The- 
oretical work on billiards where the curvature is caused by additional 
forces, both a constant force field and a parabolic potential, have been 
done (Dullin, 1998), and the stability conditions of one-periodic and 
symmetric two-periodic orbits were found. 

We investigate classical dynamics of curved-trajectory billiards that 
would have been hyperbolic if the particle had moved in straight lines. 
We introduce the curvature by applying external force fields on the parti- 
cle. We investigate two special cases: (a) a uniform force field (trajecto- 
ries are sections of parabolas), and (b) a parabolic potential (trajectories 
are sections of ellipses if the potential is attracting, and sections of hyper- 
bolas if the potential is repulsive). These perturbations can lead to the 
formation of elliptical orbits in hyperbolic billiards. We experimentally 
demonstrate these effects using atom-optics billiards. 



A periodic trajectory exists in the billiard if there is a "fixed point" 
in its Poincare map, T. In general, the map is nonlinear and cannot be 
written as a 2 x 2 matrix. In order to analyze the stability of an orbit, we 
introduce the "stability matrix" , M, which is a locally linearized version 
of T, and controls the evolution of an initially infinitesimal deviation 
from the starting point. Periodic orbits can be either stable or unstable, 
in the sense that an initial orbit, starting at (so + Ssq,pq + 5po), with 
8sq and Spo infinitesimally small, will remain close to the periodic orbit 
starting at (so,po) or become distant from it, when the map is iterated. 
After N iterations, the periodic orbit returns to its initial value, and the 
stability of the orbit is determined by the eigenvalues of L = M N , which 
are given by: 
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Periodic Orbit Stability 




(13) 
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If |trL| < 2 then both eigenvalues have absolute value equal to unity and 
are complex conjugates of the form A± = exp [±«/3]. The deviations from 
the periodic orbit oscillate around zero, but remain bounded. Hence, the 
orbit is stable or "elliptical". If |trL| > 2 then both eigenvalues are real, 
positive and reciprocals of each other, |A±| = exp [±7]. Because of the 
positive exponent deviations will grow exponentially and the orbit is 
unstable or "hyperbolic" . 

We limit ourselves to a two-bounce orbit, with impacts on both sides 
at normal incidence. In this case, the trace calculation leads to the 
well known geometrical stability conditions, in terms of p, the distance 
between the two scattering points, and i?i and R2, the corresponding 
radii of curvature, defined as positive if the particle scatters from a 
concave surface: 



Note, that these are the same condition as geometrical optics gives for 
the stability of a cavity. 

10.2 A Constant Force Field 

We first see how the presence of a uniform force field (constant ac- 
celeration g) affects the stability. The motion in between bounces is no 
longer in straight lines, but along sections of a parabola. For simplicity 
we assume a two-periodic orbit along the force field. We assume that 
the particle possesses enough kinetic energy to access the entire billiard 
and further limit ourselves to an orbit which is in the direction of the 
force field F. In this case the linearized Poincare map is given by: 



where p\ and P2 are the magnitudes of the velocities at encounters with 
the wall and t is the time it takes the particle to move from one bounce 
to the next. This trace is very similar to the trace without a force field. 
There are in general four solutions to the equation |trL| = 2 (just as 
without the force field), namely t = 0, t = ^, t = ^ and t = ^ + ^. 
Since p = —\gt 2 + pit, an interpretation identical to the one without 
the force field is also valid, provided that instead of the geometrical radii 
of curvature R, we define a pair of "effective" radii given by: 




(14) 




(15) 




i = 1,2 



(16) 
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Figure 10. Poincare surface of sections for classical trajectories of Rb atoms inside 
the half Bunimovich billiard, with the parameters used in the experiment, (a) with 
arc up: a large island of stability appears around the vertical two-periodic orbit, 
(b) with arc down: completely chaotic phase-space. Similar phase-space diagrams 
with no gravity reveal a similar completely chaotic phase-space. Insets show typical 
trajectories 
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Figure 11. (a) The fraction of remaining atoms as a function of time for Half a 
Bunimovich stadium in a constant force field (gravity). : With the arc up. ■: 
Arc down. The existence of a stable island in phase-space causes a slowing down 
in the decay. Inset: Measured light distribution of the billiard (205 fim high), with 
side sections removed, (b) Fraction of remaining atoms in a dispersing billiard as a 
function of time when a nearly-parabolic attractive potential is induced by a standing 
wave beam of variable power. ■: mW, : 10 mW, A: 20 mW and V: 30 mW. 
The formation of a stable island is seen as a slowing down in the decay. Inset: The 
billiard. 



Using the effective radius the stability can be found from Eq. 14. 
Note that these effective radii correspond to a new pair of center points, 

which are now no longer the geometrical centers but are shifted by the 

— R 2 

force field to: d = C + ^Q-^- The shifted centers correspond to points 
where orbits starting on the osculating circle, and perpendicular to it, 
will cross. 
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In the shifted center picture it is now clear that a uniform force field 
can make stable islands in a hyperbolic billiard of the Bunimovich kind. 
We demonstrate this experimentally in half a Bunimovich stadium which 
is obtained by cutting it with a straight wall along its short symmetry 
axis (see insets in Fig. 11a). This preserves the hyperbolicity of the 
billiard without a force field. As a force field we use gravity. The half 
Bunimovich stadium is placed with gravity along its symmetry axis and 
the arc up. The expression for the gravity-shifted centers indicates that, 
for atoms with low enough mechanical energy, the center of the upper arc 
is shifted to outside the billiard, and thereby the orbit becomes elliptical, 
and an island of stability appears (the center of the lower line is always 
at infinity). Alternatively if the billiard is turned upside down, i.e. the 
sign of g is changed, the gravity-shifted center is inside the billiard and 
it stays hyperbolic. We stress that although the stability of this orbit 
is determined in the shifted center picture, the dynamics of the entire 
billiard is not equivalent to the dynamics of a billiard without a force 
field constructed according to the shifted centers. Indeed, we will see 
new effects such as mixed phase-space and energy dependent stability, 
which would not occur in equivalent billiards without a force field. 

Our half Bunimovich stadium is composed of a semicircle of radius 
R\ = 158 fim connected to a straight section (i?2 = °°) by two 47 
fim straight sections. The velocity regime for which the vertical two- 
periodic orbit is stable is 10.5 v reco u < P2 <14 v reco u for the arc up 
(velocity measured at the lower straight section). This velocity interval 
is populated in the experiments due to the atomic velocity spread. To 
map phase-space we performed numerical simulations of the experiment. 
Fig. 10 shows the result of the simulations in the form of a Poincare 
surface of section. As seen in Fig. 10a (the arc up), a large stable island 
is generated by gravity around the orbit predicted to be elliptical . We 
also performed simulations without gravity, which revealed a completely 
chaotic phase-space, similar to the arc down picture in Fig. 10b. 

The experimental results are shown in Fig. 11a. The fraction of 
remaining atoms is detected both with the arc up and with the arc 
down (where the two-periodic orbit is always unstable). Large stability 
is clearly observed for the arc up in agreement with the above description. 
A slowing down of the decay is seen also with the arc down for t > 40 
ms. This we attribute to ~ 1% fraction of slow atoms that do not reach 
the upper section of the billiard, and are trapped in the lower arc by the 
principle described in (Wallis et al., 1992). This is verified by observing 
the same fraction of remaining atoms at t > 40 ms without the upper 
straight section. When the same was done with the arc up, leaving only 
the lower straight section, no stability was detected, as expected. 
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10.3 



Parabolic Potential 



We next consider what happens if the atoms move in an attractive 
harmonic potential, so the trajectories in-between scattering from the 
walls are sections of ellipses. For simplicity we view only what happen to 
orbits along the radial direction of the potential. The motion in between 
bounces follows the equation of motion x = —ojx, and together with the 
initial conditions x(t = 0) = X\ and p x (t = 0) = pi, we calculate for the 
trace: 



trL = 2 



1 + 2 



P1P2 
uj 2 RiR 2 



1 sin 2 ujt - 



1 



i?2 R\ 



sin 2ojt 



• (17) 



The stability can again be found from Eq. 14 using shifted centers or 
effective radii. The center of the coordinate system is at the center of 
the potential, and the shifted center will be at 3 : 



d = ±Ci/y/l + (u>R i /p i ) 2 , (18) 

where the minus refers to Ri < —p 2 / (u) 2 \xi\). X{ is the position of 
the scattering point. From the shifted centers we see that, in contrast 
to the uniform force field, a two-periodic orbit can become elliptical 
even if i?i and i?2 are both negative, since the R corresponding to the 
shifted centers can be positive even though the geometrical R is negative 
(Ri < —p 2 / (oj 2 \xi\)). This means that islands of stability can occur 
around two-periodic orbits in dispersing billiards and Sinai billiards. 
It is also easily seen that a potential placed asymmetric on the long 
symmetry axis in a Bunimovich stadium can change the stability of the 
orbit. 

We observed islands of stability induced by an attractive parabolic 
potential experimentally in a dispersing billiard. The billiard is 220 \xm 
high and consists of four convex arcs with the same curvature for each 
opposite pair (see inset of Fig. lib). Two of the arcs have a small 
radius of R = —130 /im and the other two are very weakly curved 
with a radius of R = —1 cm. The parabolic potential is generated 
by overlapping a standing wave with the billiard (an additional retro- 
reflected laser beam). The standing wave is red-detuned 2.5 nm from 

3 Defining the shifted center of the osculating circle as the point where orbits starting perpen- 
dicular to it and displaced an infinitesimal distance along it cross, leads to that every circle 
have two centers (orbits are ellipses so they cross the x-axis twice). The stability of the orbit 
can be determined from one of the centers, and the ± makes sure that it is the correct one 
that is used. 
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resonance and has a radial Gaussian profile with dimensions larger than 
the dimensions of the billiard (wq = 244 fim), so the potential inside the 
trap can be approximated with a parabolic potential (the effect of gravity 
is then simply to shift the center of the potential). The decay curve was 
measured for several powers of the red-detuned laser beam. The results, 
presented in Fig. lib, clearly show a slowing down of the decay as the 
power grows, indicating that an island of stability is formed and grows 
in size for stronger attractive potentials. When the hole was placed 
in the bottom of the billiard, we observed a fast decay independent 
of the strength of the parabolic potential, indicating that the island is 
indeed formed around a vertical trajectory. It was also verified that no 
atoms were trapped in the red detuned standing wave by itself for the 
powers we used. We performed numerical simulations that confirmed 
that the island is centered around the two-periodic orbit we predict 
should become elliptical. 

11. Summary 

A new experimental system, the "atom-optics billiard", was intro- 
duced in this lectures, and several non- linear dynamics experiments were 
reviewed. First, the validity of the system as an experimental realization 
of the billiard model was assessed by the observation of integrable and 
chaotic motion of cold atoms. Next, we introduced a controlled amount 
of scattering to the atomic motion by exposing the confined atoms to 
a series of resonant light pulses. The measured decay curves for an 
originally stable case show that a complete destruction of the stability 
occurs. 

Next, a numerical and experimental observation of the formation of 
islands of stability in a billiards with soft walls was presented. An island 
of stability was demonstrated in a tilted Bunimovich stadium billiard, 
around the periodic trajectory that connects the points where the big 
semicircle joins the straight lines. Our results show that the appearance 
of a stable island in a soft- wall billiard is in general accompanied by areas 
of "stickiness" surrounding it. Although the size of the stable island is a 
sensitive function of softness, not always showing a monotonic behavior, 
the "sticky" trajectories modify this behavior and, in general, introduce 
a monotonic slowing down in the decay of the number of trapped parti- 
cles, as the walls are made softer. 

Finally, we analyzed the effect of adding potentials on billiards, and 
found that they can cause elliptical orbits in otherwise hyperbolic bil- 
liards. We experimentally demonstrated these effects and also found 
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that an unstable cavity placed in gravity can serve as a velocity selective 
cavity which is stable only for a narrow velocity group. 

Our research provides a well-controlled system for investigations of 
classical and quantum chaos, and a great degree of control on the struc- 
ture of phase space. Combined with the ability to perform precision 
spectroscopy, our system may have both fundamental and practical im- 
portance: From the fundamental point of view, the spectroscopy can 
serve as a very sensitive probe for the dynamics (Andersen et al. , 2003) , 
and in particular shed some light on the quantum-classical crossover (An- 
dersen et al., 2004a; Andersen et al., 2004b). From the practical point 
of view, our system provides a playground for studying the connection 
between dynamics and decoherence, which is of outmost importance in 
the emerging field of Quantum computation. In addition, understand- 
ing this connection will help develop methods to increase the coherence 
time for trapped atoms, and thereby further improve the precision and 
accuracy of spectroscopic measurements. 
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